Second order invariants and holography 



Luca Bonanno^^, Gerardo lannone'' and Orlando Luongo-^^^ 
^Dip. di Fisica, Universita di Roma "La Sapienza", 1-00185 Roma, Italy, 
^ Institut far Theoretische Physik, J. W. Goethe - Univ. Max-von-Laue-StraQe, 1 60438 Frankfurt am Main, Germany, 
^ Dip. di Fisica, Universita di Ferrara and INFN, via Saragat 1, 44^00 Ferrara, Italy, 
^ Dip. di Scienze Fisiche, Universita di Napoli "Federico 11", Via Cinthia, Napoli, Italy, 
^Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, and 
Dip. di Fisica, Universita di Salerno and INFN, 1-84081, Salerno, Italy. 



o 



Motivated by recent works on the role of the Holographic principle in cosmology, we relate a 
class of second order Ricci invariants to the IR cutoff characterizing the holographic Dark Energy 
density. The choice of second order invariants provides an invariant way to account the problem of 
causality for the correct cosmological cutoff, since the presence of event horizons is not an a priori 
assumption. We find that these models work fairly well, by fitting the observational data, through a 
combined cosmological test with the use of SNela, BAO and CMB. This class of models is also able 
to overcome the fine-tuning and coincidence problems. Finally, to make a comparison with other 
recent models, we adopt the statistical tests AIC and BIG. 

PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es 
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I. INTRODUCTION 

Ever since the discovery of the Dark Energy (DE) [I1-[H , 
the fundamental problem of understanding the nature of 
the positive acceleration of the universe appears lacking 
of some ingredients in the framework of General Relativ- 
ity (GR). In addition, observations strongly suggest that 
matter in the universe is dominated by a non-baryonic 
(cold) dark matter (DM) which seems to evolve sepa- 
rately from the unknown force driving the acceleration 
0, [^. Despite neither DE nor DM have been directly 
detected, they represent the most relevant amount of en- 
ergy and matter in the universe, being almost the 70 and 
25 percent of the total content, respectively. A theory 
propounded to explain their effects and to describe their 
natures is strictly necessary in cosmology, although, u 
to now, it has not been still definitively formulated 
Probably, the most common paradigm trying to include 
both DE and DM, is the so-called ACDM model f^. In 
this model, a cosmological constant term, A, is included 
within the Einstein equations and the total matter den- 
sity is assumed to be the sum of baryonic and DM den- 
sities, i.e. p = Pb + Pdm- This model fits excellently and 
with high-precision all observational data, meanwhile it 
provides a remarkably small number of cosmological pa- 
rameters Q. The effects of A are hidden in the defini- 
tions of the pressure, p, and of matter density, p, i.e. 
p ^ p + P ^ P — ^g'j thus, by assuming hereafter 
a (flat) Friedman-Robertson- Walker (FRW) metrics 



a{tf [drV(f - kr^) + sin^ Odcj)'^ + dO"^ 



I.e. 

we 



can easily write down the correspondent Friedmann equa- 
tions: 
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where G is the gravitational constant, and i?o the radius 
of the universe. If p and p are known (or equivalently 
H is known), the above coupled equations describe the 
dynamics of the universe. 

In spite of its simplicity and its triumph in fitting data, 
the ACDM model suffers from several shortcomings Q, 
as e.g. the problem of fine tuning and the problem of co- 
incidence. The problem of fine tuning is connected to the 
large difference between the value of vacuum energy in 
quantum field theory and the observed value of the cos- 
mological constant . The problem of coincidence is re- 
lated to the incredibly small differences between the den- 
sities of DE and DM, which are supposed to evolve differ- 
ently as the universe expands. Moreover, one can notice 
that at large enough energies (typically of the order of 
the Planck scale, M^i = y/hc/G - 1.22 x IQi^GeV/c^), 
the ACDM model is supposed to be not reliable anymore; 
in other words, it only provides a limited description ffT| 
of the early stages of the universe. As a consequence, 
the crucial role played by inflation still remains an effec- 
tive approach without any basis in a fundamental theory. 
Nevertheless, a Quantum Gravity (QG) theory seems to 
be necessary to explain these stages (7| , having as a limit 
the ACDM model; that is clear if one believes that GR is 
a particular case of a more fundamental theory. This is 
the underlying philosophy of what are usually referred to 
as modified theories of gravity, concerning the modified 
Einstein-Hilbert (EH) action jfj, 13 • Besides fundamen- 
tal physics motivations, cosmologists acquired a huge in- 
terest in all of these theories, thanks to the possibility to 
reach a unified scheme. Among the various alternatives, 
crucial advances have been carried out in the studies of 
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Black Hole theory and String theory jf 4 |. 

A very intriguing concept, within the framework of 
GR, is the so-called holographic principle (HP), which 
provides some clues for solving the modern theoretical 
and observational problems without directly modifying 
the EH action. As it has been pointed out in QG, the 
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entropy of a system does not scale with the volume of 
that system, but with the area of its surface [l5| . Start- 
ing from the above consideration, the HP postulates that 
the maximum entropy inside a region is not extensive, 
but grows as the area of the surface. Therefore, the total 
number of independent degrees of freedom should scale 
with the surface area (in Planck units) as well. In par- 
ticular, the principle invokes that L^pvac < ^Mp^, where 
Mpi is the Planck mass and p^ac is the vacuum energy 
density of a system of size L. Actually, the HP repre- 
sents a new basic principle for both QG and GR, being 
supported by an effective quantum field theory. In order 
to include this principle in a real cosmological scenario, 
one needs to choose the correct cosmological length scale 
L, which is not a priori known. 

Since the cosmological system is the universe, the 
DE should be associated with the scale density, namely 
puE = Px oc L. Writing the fraction of DE density in 
the form Q.x = ^^'i^ jji , we can imagine different holo- 
graphic approaches by choosing different IR cutoffs. For 
instance, for the case of the ACDM, if the Hubble pa- 
rameter H is taken to be the characteristic scale of the 
universe, then it is natural to postulate A^^/^ oc H; un- 
fortunately this scenario fails in reproducing the positive 
acceleration of universe. Nevertheless, from the HP prin- 
ciple, it would be possible to solve both the problems 
of fine tuning and coincidence just by introducing the 
correct length scale. In order to fix such a scale, many 
approaches have been investigated in literature. In the 
so-called Holographic Dark Energy model (HDE) [l6l . [l7j 
the future event horizon of the universe characterizes the 
length L, i.e. 



L — a 



da' 
Ha^ 



(2) 



An alternative is the so-called agegraphic model (ADE) 
P^I - |20| , where the IR cutoff is the conformal age of the 
universe, 



da' 
a'^H 



(3) 



Unfortunately, the choice of a length scale is neither so 
easy nor arbitrary, since it leaves unclear some conceptual 
problems. As pointed out by Cai 2l[, a drawback, con- 



cerning causality, appears in the above scenarios: how is 
it possible that a local quantity, as e.g. the DE, could be 
explained by global concepts, derived from the physics of 
space-time? Is there a mechanism allowing a local quan- 
tity to be determined by a global one j2^j ? A crucial issue 
is that the above cosmological lengths are originated from 
an expanding universe, which is an a priori assumption 
and not the result of a certain model, as it should be. In 
order to solve this problem and inspired by the HP, we 
suggest to choose as length scales only those quantities 
which are invariant under geometrical transformations, 
avoiding the causality issue. This allows to solve both 
the coincidence and fine-tuning problems, as well. The 



invariants of curvature, derived from the Ricci scalar i?, 
the Weyl tensor G^^^s and the Riemann tensor R^u^s, 
are based on such theoretical implicit assumptions and 
represent geometrical invariants. Gao [235 suggested, for 
the first time, that the DE density could be proportional 
to the first order invariant 



Px oc R. 



(4) 



In this work we extend the analysis by Gao to the case 
of second order (independent) invariants. Notice that 
a second order invariant is proportional to the inverse 
fourth power of an IR cutoff, then indicating with the 
i-th independent second order invariant, our assumption 
should be formulated as follows 

Px oc ^/\^. (5) 

The purpose of the present article is to study the models 
originating from these invariants and their implications 
on the FRW spacetime. As shown below, these models 
fairly overcome the problem of the length scale, providing 
a good solution for the DE problem 171] and indicating a 
good agreement with the theoretical predictions. 

The paper is organized as follows: in Section II we 
develop the theory of the independent second order in- 
variants for a flat FRW universe and we analyze the out- 
comes of the assumption ([S|) in cosmology. In Section HI 
we flt our models with the cosmological data from Super- 
novae la (SNela), Baryonic Acoustic Oscillation (BAO) 
and Cosmological Microwave Background (CMB). In Sec- 
tion IV, we provide a comparison with other relevant cos- 
mological models and we make use of model independent 
statistical tests, i.e. the AIC and BIG. Finally in Section 
V we present the conclusions. 



II. SECOND ORDER HOLOGRAPHIC 
INVARIANTS 

The basic purpose of this work is to invoke the HP and 
to relate the second order curvature (independent) invari- 
ants to the DE density (as already discussed, the case of 
the first order curvature invariant has been extensively 
treated by Gao [H-ii]). 

This can be easily performed by writing down the ex- 
plicit form of the second order geometrical invariants 
for a FRW metric and by using eq. ([S]). As pointed 
out in Ref . [27| . among the 14 curvature scalar invariants 
[26j . the most interesting ones are the Kretschmann, the 
Chern-Pontryagin and Euler invariants. We can write, 
for every spacetime (28l-l32|. their expressions as follows 

Ki = RapysR"'^^^ , 

K2 = [*RU^sR"^''' , 

K3 = [*R*U^sR''^^', (6) 



^ Usually the model is referred to as the Ricci Dark Energy (RDE) 
model. 
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where the stars indicate the correspondent dual counter- 
parts. From the first Matte-decomposition of the Weyl 
tensor, it is easy to get [sl] 
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Therefore, Ki, K2 and K3 can be expressed as follows 
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2 , 
-R^ = 



— /l + ^RafiR^^ — T^R 



From the above equations it is straightforward to infer 
the explicit expressions of the second order invariants in 
a FRW universe: 
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respectively. In order to fix the free parameters of our 
two models, we need to find for which values of a, the 
conditions 50 < and w ^ — 1 at redshift z = are ful- 
filled. The corresponding intervals are a € [0.042, 0.052] 
and a £ [0.039, 0.047] for modi and mod^ respectively. In 
figures [H [21 [3] and [3] we plot q{z) and w{z) as functions 
of the redshift for both modi and mods. As shown in 
the pictures, if we set a = 0.050 for modi and a = 0.046 
for mods we get the expected accelerated behavior of the 
universe, being q ~ —0.6 and w ~ — 1. In particular, the 
results obtained with the two models suggest that the 
universe is accelerating since z ^ 0.4. 
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FIG. 1. Deceleration parameter g as a function of the redshift 
z for modi (in this case we use conventionally a — 0.050). 
Note that q is negative for z < 0.4, providing an accelerated 
behavior of the universe since z = 0.4; this is quite in agree- 
ment with the prediction of ACDM which refers to z ~ 0.77. 



The HP postulate reads 
3a 

Px = 



z = l,2,3; 



(10) 



where a is an dimensionless constant and should not be 
confused with the tensorial index. 

By combining together eqs. PU)) and ([IJ and by using 
the above expressions for Ii and I3 (the only two non- 
trivial invariants), we obtain two differential equations, 
each one providing the temporal evolution of the Hubble 
parameter. A first interesting step consists in solving the 
associated differential equations numerically for both Ii 
and I3, finding the correspondent acceleration parame- 
ters, by using the definition 



(11) 



and the effective barotropic parameters, expressed by 

1 d 
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(12) 



for both the models. 

From now on, we will refer to the cosmological models 
arising from the invariants Ii and I3 as modi and mod^, 




0=0.0502 



FIG. 2. The barotropic factor w = - as a function of the 
redshift z for modi (in this case we use a = 0.050). Notice 
that TO = — 1 at z = 0, behaving as a cosmological constant 
at low redshift. 



Once determined q{z) and w{z) we can characterize 
completely the kinematics of the models. In the following 
section we will check the agreement of our models with 
the observations. In particular, we will make use of three 
independent tests: SNela, BAO and CMB. 
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FIG. 3. Deceleration parameter g as a function of the redshift 
z for modz (in this case we use a = 0.046). Again q is negative 
for z < 0.4, providing an accelerated behavior of the universe 
since z — 0.4 until now. 
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FIG. 4. The barotropic factor w as a function of the redshift 
2 for mods (in this case we use a = 0.046). Notice that again 
w = —1 at z — 0, behaving in this case too as a cosmological 
constant term at low redshift. 



III. COSMOLOGICAL CONSTRAINTS 

In this section we perform an experimental combined 
procedure, by using a fairly typical combination of kine- 
matical data. In particular we employ the three most 
common fitting procedures: SNela, BAO and CMB; the 
first two concern low redshift data sets spanning from 
z = to z ^ 2, while the third is a higher redshift test, 
since it is performed at z ^ 1000. 

It is well established that standard candle data from 
SNela are indicators of distance, able to fit the corre- 
spondent luminosity distances for a particular class of 
models, by considering the distance modulus and the cor- 
respondent redshift z. The problem of the systematics, 
which usually afi^ects these measures, can be avoided or 
at least reduced, by the use of the most recent updated 
Union 2 compilation [s^J, instead of other older samples 
[ssj . Then, associating to each Supernova modulus fi the 
corresponding la error, denoted by cr^, we can perform 
directly the experimental analysis. To this purpose, let 



us rewrite the distance modulus 



— 25 + 5 logj 



di 



where dL{z) is the luminosity distance, defined by 



dL{z) =c(l + z) 



dz' 
W) 



(13) 



(14) 



Since the likelihood function C is related to the chi-square 
statistic, i.e. C oc exp(— x^/2), we constrain the free 
parameters of a model, by minimizing the quantity 



E 



the 
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(15) 



The second test that we perform is related to the ob- 
servations of large scale galaxy clusterings, which pro- 
vide the signatures of the BAO [s^. In particular we 
use the measurement of the peak of luminous red galax- 
ies observed in Sloan Digital Sky Survey (SDSS), usually 
denoted by A and defined as follows 



-HizBAoY 
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dz 



(16) 

with ZBAO = 0.35. In addition, the observed A is 
estimated to be 



Aobs = 0.469 



0.95 
OJS 



-0.35 



(17) 



with an error aA = 0.017. In the case of the BAO mea- 
surement we should perform the following minimization: 



y^BAO 



A 



<JA 



(18) 



Finally, concerning the CMB, we first analyze the so- 
called CMB shift parameter 



R 



Ho 
H{z) 



dz . 



(19) 



Since this standard ruler is fixed by the sound horizon 
at decoupling (^doc = 1091.36 |37|]), it gives a complemen- 
tary bound to the SNela data and BAO as well. For the 
R parameter, the observed value is Rq^ = 1.726 ±0.018, 
as inferred from the WMAP 7 data [3J|. Hence, the cor- 
respondent CMB constraints are given by minimizing the 
chi square 



XcMB — 



R — R, 



obs 



(20) 



Moreover, we notice that, differently from SNela, BAO 
and CMB do not depend on Hq. 

Finally, by combining the minimization procedures for 
SNela, CMB and BAO, we constrain the parameters of 
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QmiSN) Q.^{BAO) QmjCMB) 

0.240 ± 0.070 0.234 ± 0.014 0.266 ± 0.064 

ajSN) a{BAO) ajCMB) 

0.046 ±0.016 0.044 ± 0.010 0.025 ± 0.016 



TABLE I. Summary of the results for modi; in this case 
we have XsNeia = 1-021, Xbao = 1-001, Xcmb = 1-001, 

i^^,(mean) = 0-247 ± 0.070, Qf^ean) = 0.038 ± 0.016. 



Vlm{SN) Q.m(BAO) Q.m(CMB) 

0.260 ± 0.090 0.312 ± 0.080 0.343 ± 0.060 

a{SN) a{BAO) a{CMB) 

0.042 ± 0.010 0.042 ± 0.002 0.024 ± 0.0015 



TABLE IL Summary of the resuhs for mods. The reduced chi 
squared are x%Neia = 1-021, Xbao = 1-001, Xcmb = 1-003; 

the mean values flm,{rnean) = 0.305±0.090, a^rnean) ~ 0.036± 

0.015. 



our two models. The numerical results are summarized 
in the following tables: 

The results show that the theoretical predictions dis- 
cussed in Sec|TT] are in agreement with the experimental 
ones. The mean values of a in the tables are indeed 
included in the intervals of values found in Sec|lll The 
CMB measurements are usually smaller if compared with 
the SNela and BAO ones. Notice finally that for the 
SNela we had to fix a value for Ho (Hq ^ 2.33-^^ s~^), 
while for CMB and BAO it is not necessary. 



IV. COMPARISON WITH OTHER MODELS 

In the Introduction we pointed out that the ACDM 
paradigm appears to be the favorite fitting model among 
a large number of possibilities. This should depend on its 
small number of parameters. In particular, for a flat cos- 
mology, the only parameter involved is the mass density 

- 

Therefore, one can ask if there is a real necessity to go 
beyond this approach, by considering other frameworks. 
This question suggests the requirement to find a test able 
to compare different cosmological models, in order to se- 
lect the "best" one. At the same time, such a test should 
also be model independent. 

A good choice is represented by the so-called Akaike In- 
formation Criterium (AIC) and BIC ^] tests, which 
are two of the most model-independent statistical meth- 
ods for comparing different models. Moreover, since their 
first use by Liddle [40] , these tests became a standard di- 
agnostic tool [4ll443{ of regression models (43 - l47| . 

The idea behind AIC is based on postulating two dis- 
tribution functions, namely f{x) and g{x\0): f{x) is as- 
sumed to be the exact one, while g{x\6) approximates the 
former through a set of parameters denoted by 6. Hence, 



once given f{x) and g{x\0), there exists only a set of Om.in. 
minimizing the difference between g{x,9) and f{x) [48'|. 

However, without going into details, we only note that 
the AIC value for a single model is meaningless since 
the exact model function f{x) is unknown. Therefore, 
the quantities of interest are the differences AAIC = 
AIC — AlCmin, calculated over the whole set of models. 
The generic AIC is given by 

AIC = -21n/:,„ai; ± 2k. (21) 

A very similar criterion was derived by Schwarz psj in a 
Bayesian context (see [i^ and references therein). 

The BIC test provides 

BlC^-2\n£raax + k\nN, (22) 

where, for both AIC and BIC, Cmax is the maximum 
likelihood, k is the number of parameters, and TV is the 
number of data points used in the fit. If the errors are 
Gaussian, then Xmin = ^2lnC„iax and the difference in 
BIC can be simplified to ABIC = Ax^„,„ ± Afc In N. 

We adopted both the AIC and BIC tests for different 
models. 

A first natural extension of the ACDM paradigm is the 
wCDM model, also called quintessence model; it provides 
an EoS of the form p = wp with a negative barotropic 
factor, whose origin is related to the coupling of the Ricci 
scalar with a not evolving scalar field (j). Since this model 
is strongly dependent onw, it is called wCDM in analogy 
with the ACDM. 

The theoretical explanation about the origin of the 
scalar field giving a negative EoS remains unsolved in 
the wCDM, then a varying quintessence has been pro- 
posed: if w evolves with the redshift z, i.e. w = w[z], 
the scalar field evolves as well. So the origin of the latter 
should be found thermodynamically or in other ways. 

One of the more recent and intriguing varying 
quintessence model has been proposed by the Chevallier, 
Polarsky, Linder (CPL) [i^. This parametrization sug- 
gests that w[a] = wo+wi{l — a). Assuming a = (l+z)~^, 
one gets w[z] = wq + wij^, which for low and very high 
redshifts becomes constant, i.e. w(z — > 0) = wo and 
w[z — >■ oo] = Wo + Wi. 

We perform the AIC and BIC tests for our two models 
{modi and mod^) and for the models discussed above. 
Moreover, we also include the Ricci DE model (RDE), 
studied by Gao [1^. The normalized Hubble rates, E = 
■jj- , for these models are written below (notice that there 
is not an analytic expression of E for the models modi 
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and modj): 

EaCDM = \/^rniz) + 1 — flm , 

E^CDM = ^a«(2) + (l-r!„0(l + z)3(i+-), 

(23) 

EcpL = v/a„(2) + (i-ri„0/(z), 



Eb.de = \U — n^iz) + nf{i + zf-l , 

V z — a 

where /(z) = (1 + cxp |-^| and 

The results of the tests are summarized in the table 
below, in which we report, for each model, the number 
and the names of the free parameters, the s-nd the 
values of BIG and AIC: 



Model 


N. Par. (k-1) 


Par am. 


'X.min 


AB/C 




kCDM 


1 




557.40 








Quint. 


2 




557.28 


6.20 


1.88 


CPL 


3 


Wo, Wa 


557.45 


12.69 


4.05 


RDE 


2 




573.72 


22.64 


18.32 


modi 


2 




568.86 


17.79 


13.46 


mod-i 


2 


Vim, a 


568.91 


17.84 


13.51 



where we have used: Q,^ = 0.254 ± 0.038 for ACDM; 
= 0.325 ± 0.049 and w = -1.18 ± 0.18 for wCDM; 
a™ = 0.246±0.034,^i(;o = -0.91±0.12 and Wa = -0.32± 
0.04 for CPL *53,l5ll; = 0.310 ±0.052 and a = 0.380 ± 
0.049 for the Ricci DE. 

The results have been obtained through a direct anal- 
ysis of each model, following the combined procedure ex- 
plained in Section III. We can conclude that our models 
are disfavored by the AIC and BIC analysis if compared 
with ACDM and wCDM. It is clear that ACDM remains 
the favorite model. On the other hand, the results ob- 
tained for the Ricci DE are even worst. This seems to 
suggest that using higher order invariants is better than 
using invariants of lower order, encouraging further inves- 
tigations in our direction. Concerning the CPL model, 



it is interesting to notice that, to a slight variation of 
Wo and Wa, it corresponds a large variation of AIC and 
BIC. This is in agreement with the fact that CPL is a 
three parameters approach, then it appears disfavored if 
compared with modi and mods. 

V. CONCLUSION 

In this work we have investigated the possibility to re- 
late the DE term with the second order independent in- 
variants li . This has been performed in analogy with the 
work of Gao (23j , who proposed a first order invariant ap- 
proach, namely the Ricci DE. The validity of our idea is 
based on the holographic principle, which requires the ex- 
istence of a cut-off scale, characterizing the dark energy. 
Moreover, the latter must scale with the inverse square 
of the cut-off (which means px oc \/Ti). After solving 
numerically the Friedmann equations for the two non- 
trivial invariants in a flat FRW universe, we tested the 
two resulting models, namely modi and mod^. Hence, we 
first analyzed the kinematics of the two models in terms 
of the acceleration parameter q{z) and in terms of the 
barotropic factor w{z), showing that they are compati- 
ble with the conditions q < and latz = 0. 
Then, we developed a test combining SNela, BAO and 
CMB in order to fix the cosmological parameters; again 
in this case we found a good agreement with the exper- 
imental results. It is also important to notice that the 
two models, arising from different invariants, show an 
incredibly similar behavior. 

A robust evidence of the validity of the models is given 
by the AIC and BIC tests, where we get results bet- 
ter than those obtained with the Ricci DE approach. 
Therefore, our models are compatible with the holo- 
graphic principle and seem to behave better than the 
first order invariant proposed by Gao. This encourages 
to study higher order invariants, what will be done in 
future works. 
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